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Abstract. We give some explicit calculations for stable distributions 
and convergence to them, mainly based on less explicit results in Feller 
[3]. The main purpose is to provide ourselves with easy reference to 
explicit formulas and examples. (There are probably no new results.) 



1. Introduction 

We give some explicit calculations for stable distributions and convergence 
to them, mainly based on less explicit results in Feller [3]. The main purpose 
is to provide ourselves with easy reference to explicit formulas and examples, 
for example for use in [8]. (There are probably no new results.) These notes 
may be extended later with more examples. 

2. Infinitely divisible distributions 

We begin with the more general concept of infinitely divisible distribu- 
tions. 



Definition 2.1. The distribution of a random variable X is infinitely divis- 
ible 
that 



ible if for each n ^ 1 there exists i.i.d. random variable Y^"'\ . . . , Y^"^ such 



x = y/"^ + ...yi"). (2.1) 

The characteristic function of an infinitely divisible distribution may be 
expressed in a canonical form, sometimes called the Levy-Khinchin repre- 
sentation. We give several equivalent versions in the following theorem. 

Theorem 2.2. Let h{x) be a fixed bounded measurable real-valued function 
on M such that h{x) = x+0(x^) as x — )• 0. Then the following are equivalent. 

(i) ip{t) is the characteristic function of an infinitely divisible distribution. 

(ii) There exist a measure M on M such that 

/oo 
(l A |xr^) dM(x) < cx) (2.2) 
-oo 

and a real constant b such that 

/ , e'*"" - l-^th(x) ,,,, ,\ , , 

ip{t) = exp(ibt+ 2 —dM{x)], (2.3) 

^ ./ — ^ X / 
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where the integrand is interpreted as — at x = 0. 
(iii) There exist a measure A on R \ {0} such that 



/oo 

A 1) dA(a;) < oo (2.4) 
-GO 

and real constants a ^ and b such that 

ip{t) = exp(ibt - iat^ + f (e'*^ - 1 - ith{x)) dA(x)) . (2.5) 

(iv) There exist a bounded measure K onM. and a real constant b such that 

^{t) = exp(i6t + J^^ (e^*- - 1 - ^) 1^ dK{x)) , (2.6) 

where the integrand is interpreted as —t^/2 at x = 0. 
The measures and constants are determined uniquely by Lp. 

Feller [3, Chapter XVII] uses h{x) = sin a;. Kallenberg [9, Corollary 15.8] 
uses h{x) = xl{\x\ ^ 1}. 

Feller [3, Chapter XVII. 2] calls the measure M in (ii) the canonical mea- 
sure. The measure A in (iii) is known as the Levy measure. The parameters 
a, b and A are together called the characteristics of the distribution. We 
denote the distribution with characteristic function (2.5) (for a given h) by 
ID{a,b, A). 

Remark 2.3. Different choices of h{x) yield the same measures M and A 
in (ii) and (iii) but different constants b; changing h to h corresponds to 
changing b to 

b:=b+ r ^(^) ~^(^) dM(x) =b+ r {h{x) - h{x)) dA{x). (2.7) 

J —oo J —oo 

We see also that 6 is the same in (ii) and (iii) (with the same h), and that 
(see the proof below) b in (iv) equals b in (ii) and (iii) when x = x/{\ + a;^). 

Proof, (i) (ii): This is shown in Feller [3, Theorem XVIL2.1] for the 
choice h{x) = sinx. As remarked above, (2.3) for some h is equivalent to 
(2.3) for any other h, changing b by (2.7). 

(ii) (iii): Given M in (ii) we let a := M{0} and dA(x) := x''^ dM(x), 
X 7^ 0. Conversely, given a and A as in (iii) we define 

dM(x) = a5o + x^ dK{x). (2.8) 

The equivalence between (2.3) and (2.5) then is obvious. □ 

(ii) <^=^ (iv): Choose h{x) = x/(l + x^) and define 

dK{x) := dM(x); (2.9) 

conversely, dM{x) = {1 + x"^) dK {x) . Then (2.3) is equivalent to (2.6). 
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Remark 2.4. At least (iii) extends directly to infinitely divisible random 
vectors in W^. Moreover, there is a one-to-one correspondence with Levy 
processes, i.e., stochastic processes Xt on [0, oo) with stationary independent 
increments and Xq = 0, given by (in the one-dimensional case) 

Ee'"^* = ip{uy = ex.p(t(^ibu-^au'^+ J (e'"^-l-in/t(a;)) dA(x))) (2.10) 

for t ^ and u G M. See Bertoin [1] and Kallenberg [9, Corollary 15.8]. 

Example 2.5. The normal distribution N{ii,a^) has A = and a = 
thus M = K = a^So; further, 6 = /x for any h. Thus, N{fi, a^) = lD{a^, fi, 0). 

Example 2.6. The Poisson distribution Po(A) has M = A = X6i and 

K = ^6i; further b = Xh{l). (Thus b = A/2 in (iv).) 

Example 2.7. The Gamma distribution Gamma(a) with density function 
x'^~^e~^ /T{a), X > 0, has the characteristic function ip(t) = (1 — it)~". It 
is infinitely divisible with 

dM(x) = axe-^, x > 0, (2.11) 

dA(x) = ax~^e~^, x > 0, (2.12) 

see Feller [3, Example XVII.3.d]. 

Remark 2.8. If Xi and X2 are independent infinitely divisible random 

variables with parameters (04,61, Ai) and (02)^2) A2), then Xi + X2 is in- 
finitely divisible with parameters (ai -f 02, 61 -|- 62, Ai -|- A2). In particular, 
if X-ID(a,6,A), then 



X = Xi + Y with Xi ~ ID(0, 0, A), Y - lD(a, b, 0) = N{b, a), 

(2.13) 

and Xi and Y independent. Moreover, for any finite partition M = IJ^j, 
we can split X as a sum of independent infinitely divisible random variables 
Xi with the Levy measure of Xi having supports in A^. 

Example 2.9 (integral of Poisson process). Let S be a Poisson process on 
M \ {0} with intensity A, where A is a measure with 

/oo 
(|x| A 1) dA(x) < 00. (2.14) 
-00 

Let X := J a;dH(a;); if we regard S as a (finite or countable) set (or possibly 
multiset) of points {^i}, this means that X := X^j^j- (The sum converges 
absolutely a.s., so X is well-defined a.s.; in fact, the sum J2\^i\>i^i 
finite, and the sum X^|^,|^i has finite expectation J^^ |x|dA(x).) Then 
X has characteristic function 

^{t) = exp( r (e'*^ - 1) dA(x)) . (2.15) 

J —00 

(See, for example, the corresponding formula for the Laplace transform in 
Kallenberg [9, Lemma 12.2], from which (2.15) easily follows.) Hence, (2.5) 
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holds with Levy measure A, a = and b = J^^h{x) dA{x). (When (2.5) 
holds, we can take h{x) = 0, a choice not allowed in general. Note that 
(2.15) is the same as (2.5) with h = 0, a = and b = 0.) 

By adding an independent normal variable N{b,a), we can obtain any 
infinitely divisible distribution with a Levy measure satisfying (2.5); see 
Example 2.5 and Remark 2.8. 

Example 2.10 (compensated integral of Poisson process). Let S be a Pois- 
son process on R \ {0} with intensity A, where A is a measure with 

/OD 
{\xf A \x\) dA{x) < oo. (2.16) 
-oo 

Suppose first that \x\ dA(x) < oo. Let X be as in Example 2.9. Then 
X has finite expectation KX = xdA. Define 

/oo 
x{dE{x)- dA{x)). (2.17) 
-oo 

Then, by (2.15), X has characteristic function 

ip{t) = exp( r (e"^ - 1 - itx) dA{x)^ . (2.18) 



Now suppose that A is any measure satisfying (2.16). Then the integral 
in (2.18) converges; moreover, by considering the truncated measures A.„ : = 
l{|x| > n~^}A and taking the limit as n — t- oo, it follows that there exists 
a random variable X with characteristic function (2.18). Hence, (2.5) holds 
with Levy measure A, a = and b = J^^{h{x) — x) dA(x). (When (2.16) 
holds, wc can take h{x) = x, a choice not allowed in general. Note that 
(2.18) is the same as (2.5) with h{x) = x, a = and b = 0.) 

By adding an independent normal variable N{b,a), we can obtain any 
infinitely divisible distribution with a Levy measure satisfying (2.16); see 
Example 2.5 and Remark 2.8. 

Remark 2.11. Any infinitely divisible distribution can be obtained by tak- 
ing a sum Xi + X2 + Y of independent random variables with Xi as in 
Example 2.9, X2 as in Example 2.10 and Y normal. For example, we can 
take the Levy measures of Xi and X2 as the restrictions of the Levy measure 
to {x : \x\ > 1} and {x : \x\ ^ 1}, respectively. 

Theorem 2.12. If X is an infinitely divisible random variable with char- 
acteristic function given by (2.5) and t E'R, then 

Ee*^ = ex.p(bt + \at^ + / (e*"^ - 1 - th{x)) dA{x)\ ^ 00. (2.19) 
^ J —00 ' 



STABLE DISTRIBUTIONS 



5 



In particular, 

Ee*^<oo^ / (e*^ - 1 - dA(x) < oo 

dA{x) < oo, t>0, 
dA{x) < oo, t <0. 




(2.20) 



Proof. The choice of h (satisfying the conditions of Theorem 2.2) does not 
matter, because of (2.7); we may thus assume h{x) = xl{\x\ ^ 1}. We 
further assume t > 0. (The case t < is similar and the case i = is 
trivial.) 

Denote the right-hand side of (2.19) by FA{t). We study several different 
cases. 

(i) . If suppA is bounded, then the integral in (2.19) converges for all 
complex t and defines an entire function. Thus F\{t) is entire and (2.5) 
shows that Ee'*^ = Fa (it). It follows that E|e*^l < oo and Ee*^ = FA{t) 
for any complex t, see e.g. Marcinkiewicz [10]. 

(ii) . If suppA C [l,oo), let A„ be the restriction Aj^^^^j of the measure 

A to [l,n]. By the construction in Example 2.9, we can construct random 
variables X„ ~ ID(0,0, A„) such that X„ X ~ ID(0, 0,A) as n ^ oo. 
Case (i) apphes to each A„, and (2.19) follows for X, and t > 0, by monotone 

convergence. 

(in). If suppA C (—00,1], let A„ be the restriction Air _ Similarly 
to (ii) we can construct random variables X„ ~ ID(0, 0,A„) with X„ ^ 
such that Xn \ X ~ ID(0, 0, A) as n ^ oo. Case (i) applies to each A„, 
and (2.19) follows for X; this time by monotone convergence. 

(iv). The general case follows by (i)-(iii) and a decomposition as in Re- 
mark 2.8. □ 



3. Stable distributions 

Definition 3.1. The distribution of a (non-dcgencratc) random variable X 
is stable if there exist constants a„ > and bn such that, for any n ^ 1, if 
Xi,X2,. . . are i.i.d. copies of X and Sn := Y17=i -^ii ^^^^ 

Sn = OnX + hn- (3.1) 

The distribution is strictly stable if 6„ = 0. 

(Many authors, e.g. Kallenberg [9], say weakly stable for our stable.) 
We say that the random variable X is (strictly) stable if its distribution 

is. 

The norming constants a„ in (3.1) arc necessarily of the form a„ = n^^" 
for some a € (0,2], sec Feller [3, Theorem VI. 1.1]; a is called the index [4], 
[9] or characteristic exponent [3] of the distribution. We also say that a 
distribution (or random variable) is a-stable if it is stable with index a. 
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The case a = 2 is simple: X is 2-stable if and only if it is normal. For 
a < 2, there is a simple characterisation in terms of the Levy-Khinchin 
representation of infinitely divisible distributions. 

Theorem 3.2. (i) A distribution is 2-stable if and only if it is normal 
N{iJ.,a'^). (This is an infinitely divisible distribution with M = a^6o, 
see Example 2.5.) 

(ii) Let < a < 2. A distribution is a-stable if and only if it is infinitely 
divisible with canonical measure 

dMjx) _ /c+x^-", x>0, ^^2) 



dx lc_|a;|^~", x < 0; 
equivalently, the Levy measure is given by 

dA(x) fc+x-"-\ x>0, 



(3.3) 



dx [c_|x|~"~^, X < 0, 
and a = 0. Here c-,c+ ^0 and we assume that not both are 0. 
Proof See Feller [3, Section XVII.5] or Kallenberg [9, Proposition 15.9]. □ 
Note that (3.2) is equivalent to 

M[xi,X2] =C+x2-" + C_|xi|2-° (3.4) 
for any interval with xi ^ ^ X2, with 

= (3.5) 

Theorem 3.3. Let < a ^ 2. 

(i) A distribution is a-stable if and only if it has a characteristic func- 
tion 



exp -7"|t|" 1 -i/3tan^sgn(t) +i6t], a / 1, 
(f{t) = \ } / \ \ (3.6) 

|^exp(^-7|t|(^l + i/3f sgn(t)log|t|j +i5tj, a = 1, 

where — 1 ^ /? ^ 1, 7 > and —00 < 5 < 00. Furthermore, an 
a-stable distribution exists for any such 0,^,^,6. (If a = 2, then (5 
is irrelevant and usually taken as Q.) 
(ii) If X has the characteristic function (3.6), then, for any n ^ 1, (3.1) 
takes the explicit form 



d J + (n - nV"),5, a^l, 

+ ^/37nlogn, a = 1. 



In particular, 

X is strictly stable <^=^ < ' ^ ' (3-8) 

1/3 = 0, a = 1. 
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(iii) An a-stable distribution with canonical measure M satisfying (3.4) 
has 

„ f(C+ + C_)£^f^ cos ^, a 7^1, 

a _ l\ -r ^a(l-a) 2 ' / ' g\ 

\(C+ + C_)f, a = l, 

(iv) // < a < 2, i/ien an a-stable distribution with Levy measure A 
satisfying (3.3) /las 

f(c+ + C-)(-r(-a)cos^), a /I, ^^.^^^ 
\(c++c_)f, a = l, 

We use the notation 80(7, 5) for the distribution with characteristic 
function (3.6). 

Proof FeUer [3, XVII.(3.18)-(3.19) and Theorem XVII.5.1(ii)] gives, in our 
notation, for a stable distribution satisfying (3.4), the characteristic function 

exp(^-(C, + C_)^i^(cos- -.sgn(t)^^sm-)ir + 
if a 7^ 1 and 

exp(^-(C+ + C_) (I + isgn(t) ^+~^- log |i|) \t\ + i6t 

if a = 1. This is (3.6) with (3.9)-(3.10) and b = b. This proves (i) and (iii), 
and (iv) follows from (iii) by (3.5). 

Finally, (ii) follows directly from (3.6). □ 

Remark 3.4. If 1 < a ^ 2, then b in (3.6) equals the mean EX. In 
particular, for a > 1, a stable distribution is strictly stable if and only if its 
expectation vanishes. 

Remark 3.5. If X - Sa(l, /3, 0), then, for 7 > and 5 e M, 

>S„(7,/3,«-f/37log7), a=l. 



Thus, 7 is a scale parameter and b a location parameter; /3 is a skewness 
parameter, and a and /3 together determine the shape of the distribution. 

Remark 3.6. More generally, \iX ^ 80,(7, ^1 then, for a > and c? G M, 



aX + ci~<l^«("^'^'"'^ + ^)' (3.14) 
' 8a(a7,^,a(5 + d- ^/37aloga), a = 1. 



8 SVANTE JANSON 

Remark 3.7. If X ~ Sa{'J,P,S), then —X ~ 80(7, — /3, —S). In particular, 
X has a symmetric stable distribution if and only if X ~ Sa(7, 0, 0) for some 
a e (0, 2] and 7 > 0. 

We may simphfy expressions hke (3.6) by considering only t ^ (or 
t > 0); this is sufficient because of the general formula 



<p{-t) = <p{t) (3.15) 

for any characteristic function. We use this in our next statement, which is 
an immediate consequence of Theorem 3.3. 

Corollary 3.8. Let < a ^ 2. A distribution is strictly stable if and only 

if it has a characteristic function 

(p{t) =exp{-{K-iX)t"), t^O, (3.16) 

where k > and |A| ^ Ktan^; furthermore, a strictly stable distribution 
exists for any such k and A. (For a = 1, tan^ = 00, so any real A is 
possible. For a = 2, tan ^ = 0, so necessarily A = 0.) 

The distribution Sa(7,,0, 0) (a ^ 1) or 81(7, 0,(5) (a = 1) satisfies (3.16) 
with 

a J X j/3Ktan^, a / 1, ,^ 

K = 7°= and A = < ^. ^ (3.17) 

Id, a = 1. 

Conversely, if (3.16) holds, then the distribution is 

|8„(7,/3,0) TOi/j7 = «V", /3 = Acot^, a/1 ^^^^^ 
I 8i(k, 0, A), a = 1. 

□ 

Remark 3.9. For a strictly stable random variable, another way to write 
the characteristic function (3.6) or (3.16) is 

<f{t) = exp(-ae'^Sn(t)7r7/2|^|a^ ^ -^q^ 

with a > and 7 real. A comparison with (3.6) and (3.17) shows that 

acos^ = «; = 7", (3.20) 

a ^ 1, , , 

^ ' (3.21) 




If < a < 1, we have < tan ^ < 00 and I7I ^ a, while if 1 < a < 2, then 

tan^ < and tan ^ = /3tan^^^^ with < 7r(2 - a)/2 < 7r/2; hence 
I7I ^ 2 — a. Finally, for a = 1, we have I7I < 1. (These ranges for 7 are 
both necessary and sufficient. For a = 1, 7 = ±1 is possible in (3.19), but 
yields a degenerate distribution X = —7a.) 
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For a 7^ 1, note the special cases (3 = 7 = and 

/3 = 1^7=|""' (3.22) 
|2-a, l<a<2. ^ ^ 

Remark 3.10. For a = 1, the general 1-stable characteristic function (3.6) 
may be written, similarly to (3.16), 

(p{t) =e^p{-{K -iX)t-ibt log t), t>0, (3.23) 

where K = -f, X = S and b = ^^j. (Thus, \b\ ^ 2K/7r.) 

Definition 3.11. A stable distribution is spectrally positive if its Levy mea- 
sure is concentrated on (0, oo), i.e., 

dA{x) = cx^""^ d.x, X > 0, (3.24) 

for some c > and a G (0,2). By (3.3) and (3.12), this is equivalent to 
c_ = and to ^ = 1, see also (3.22). A strictly stable distribution with 
characteristic function (3.16) is spectrally positive if and only if a 7^ 1 and 
A = Ktan^. 

Theorem 3.12. Let < a < 2. An a-stable random variable X ~ 
Sq;(7,/3,5) has finite Laplace transform Ee~*^ for t ^ if and only if it 
is spectrally positive, i.e., if (3 = 1, and then 



Ee-*^ =\ S^""'- J ^ (3.25) 
[expi^-jtlogt-Stj, a = l, 

Moreover, then (3.25) holds for every complex t with Ret ^ 0. 

Proof. The condition for finiteness follows by Theorem 2.12 and (3.3), to- 
gether with Definition 3.11. When this holds, the right-hand side of (3.25) 
is a continuous function of t in the closed right half-plane Ret^O, which is 
analytic in the open half-plane Rot > 0. The same is true for the left-hand 
side by Theorem 2.12, and the two functions are equal on the imaginary axis 
t = is, s G M by (3.6) and a simple calculation. By uniqueness of analytic 
continuation, (3.25) holds for every complex t with Ret ^ 0. □ 

Theorem 3.13. An stable random variable X 

^ 8(^(75/^5'^) ^s posztxve, I.e. 
X > a.s., if and only z/0<q<1, /3 = 1 and 5^0. Equivalently, X > 
a.s. if and only if X = Y + 6 where Y is spectrally positive strictly a-stable 
with < a < 1 and 6^0. 

Proof. X > a.s. if and only if the Laplace transform Ee"*"'^ is finite for 

all t ^ and Ee"*"'^ — ?■ O as t — ?■ 00. Suppose that this holds. We cannot 
have a = 2, since then X would be normal and therefore not positive; thus 
Theorem 3.12 applies and shows that ^0 = 1. Moreover, (3.25) holds. If 
l<Q;<2ora = l, then the right-hand side of (3.25) tends to infinity as 
t — > 00, which is a contradiction; hence < a < 1, and then (3.25) again 
shows that 6^0. 

The converse is immediate from (3.25). □ 
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Example 3.14. If < q < 1 and d > 0, then X ~ 80(7, 1,0) with 7 := 
(dcos 2^^!/° jg positive strictly stable random variable with the Laplace 
transform 

Ee"*^ =exp(-di"). (3.26) 

Example 3.15. If 1 < a < 2 and d > 0, then X ~ S„(7, 1,0) with 

7 := ((i| cos ^1)^^" is a strictly stable random variable with the Laplace 
transform 

Ee-*^ = exp(dr). (3.27) 

Note that in this case cos ^ < 0. Note also that Ee"*''^ — >■ 00 as t — >■ 00, 
which shows that F{X < 0) > 0. 

3.1. Other parametrisations. Our notation 80,(7,/?, (5) is in accordance 
with e.g. 8amorodnitsky and Taqqu [12, Definition 1.1.6 and page 9]. (Al- 
though they use the letters 8q,(c7, /3, 12).) Nolan [11] uses the notation S{a, j3, 7, 6 
he also defines S{a, j3, 7, do;0) := S{a, (3, 7, Si; 1) where 



So-p-ftan^, a 7^1, 
'^o-f^7log7, a = 1. 



5i := ' ''''' 2 ' " ^ (3.28) 



(Note that our S = Si.) This parametrisation has the advantage that the 
distribution S'(a, /3, 7, 5o; 0) is a continuous function of all four parameters. 
Note also that 5(a, 0,7, 5;0) = S{a,0,'y,S;l), and that when 7 = 0, (3.13) 
becomes 7X + 5 ~ 80,(7, ^, 5; 0). Cf. the related parametrisation in [12, 
Remark 1.1.4], which uses 



Ui+/37"tanf^ = <5o + /3(7"-7)tan^, a 7^ 1, ^9) 
\Si = So- ^l3j\ogj, a = 1; 

again the distribution is a continuous function of (a, /3, 7, /ii). 

Zolotarev [13] uses three different parametrisations, with parameters de- 
noted (a, /3a;, 7a:, A^,), where x G {A, B, M}; these are defined by writing the 
characteristic function (3.6) as 

^{t) = exp(AA(it7A - \tr + itcoA{t,a,PA))) (3.30) 

= exp(AM(it7M - + itcjM{t,a,pM))) (3.31) 

= exp(AB(it7B - \t\''coB{t,a,PB))), (3.32) 



where 



\t\°'-^ptan^, a 7^1, 
-;3flog|t|, a = l; 



(|t|"-i - l)^tan^, 07^1, 
-/3flog|t|, a = l; 



-B(*,a,/3):=|rp(r;f^;^("^7*)' ^^j' 

I f -Fi/31og|i|sgnt, a = l, 



(3.35) 
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with K{a) := a — 1 + sgn(l — a), i.e., 

K{a) :=h ° < " < (3.36) 

^ ^ [a-2, l<a<2. ^ ' 

Here a is the same in all parametrisations and, with /3,7, 5 is as in (3.6), 





= /3m = /3, 




(3.37) 


lA 


= Vt", 




(3.38) 


IM 




[7^ + ^tan^, a /I, 
[7A, a = 1, 


(3.39) 


\a 


= Am = 7", 




(3.40) 



and, for a / 1, 

7rivr(Q;)\ TTQ; ttq; 



tan^/3B^^j =^Atan— =;3tan — , (3.41) 
7b = 7^cos(/3b^^^), (3.42) 
Ab = Aa/cos(^b^^^), (3.43) 



while for a 



Pb = I3a = /3, (3.44) 

n tt5 

IB = 2 7^ = (3.45) 
2 2-^ 

Ab = -Aa = — . (3.46) 

TT TT 

Note that, for any a. Pb = P = and Pb = il P = il, and 

that the mapping P = Pa '-^ Pb is an increasing homeomorphism of [—1, 1] 
onto itself. 

In the strictly stable case, Zolotarev [13] also uses 

(^(t) =exp(-Ace-^'^s°W™'^/2|^|aj^ ^347^ 

which is the same as (3.19) with 

Ac = a (3.48) 
= (3.49) 

thus 1^1 ^ 1 for a ^ 1 and \d\ ^ 2/a — 1 for a > 1. We also have 

{a K(a) I -, 

, . ' (3.50) 

- arctan(27B/7r), a = 1. 

\ J^-S' 

[AB(7r74 + 7|) ' , a = \. 
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He also uses the parameters a, Ac where p := (1 + 9)/2. 

Zolotarev [13] uses Y{a, I3x-,^x-,)^x) as a notation for a random variable 
with the characteristic function (3.30)-(3.32). The distribution is a contin- 
uous function of the parameters (a, ^5,75, As). 



A stable distribution has by (3.6) a characteristic function that decreases 
rapidly as t — >■ ±00, and thus the distribution has a density that is infinitely 
differentiable. 

In the case a < 1 and /3 = 1, 8^(7, /3,(^) has support [5, 00) and in the 
case a <1 and /3 = —1, S„(7,/3,5) has support (—00, 5]; in all other cases 
the support is the entire real line. Moreover, the density function is strictly 
positive in the interior of the support, sc Zolotarev [13, Remark 2.2.4]. 

Feller [3, Section XVII. 6] lets, for a 7^ 1, p{x;a,^) denote the density 
of the stable distribution with characteristic function (3.19) with a = 1. 
A stable random variable with the characteristic function (3.19) thus has 
the density function a"^/"jo(a^^/".T; a, 7). The density of a random variable 

~ Sq,(7, 5) with a / 1 is thus given by 



with a and 7 given by (3.20)-(3.21). 

Feller [3, Lemma XVII. 6.1] gives the following series expansions. 

Theorem 4.1. (i) If x > and < a < 1, then 



4. Stable densities 



a 



VV(a-V"(x-5);a,7), 



(4.1) 




(4.2) 



k=l 



(ii) If X > and 1 < a < 2, then 




(4.3) 



fc=i 



For X < we use p{—x; a, 7) = p{x; a, —7) . 



□ 




(4.4) 



p{0; a, 2 - a) = -r(l + 1/a) sin ^ — ^ = -r(l + l/a) sin - 




r(l + l/a) 1 1 



(4.5) 



r(l/a)r(l - l/a) ar(l-l/a) |r(-l/a)|' 
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For 1 < a < 2, the distribution 80,(7, 1,0) t^^us has, by (4.1) and (3.20), 
the density at a: = 

a-V>(0;a,2-«) = " = ^-i|cos ^|^/"|r(-l/a)|-^ (4.6) 
|r(-l/a)| 2 

5. Domains of attraction 

Definition 5.1. A random variable X belongs to the domain of attraction 
of a stable distribution JC if there exist constants o„ > and 6„ such that 



d 



C (5.1) 

as n ^ 00, where Sn '■= 5^"=i Xi is a sum of n i.i.d. copies of X. 

We will in the sequel always use the notation Sn in the sense above (as 
we already have done in Section 3). All unspecified limits are as n — > 00. 

Theorem 5.2. Let < a ^ 2. A (non-degenerate) random variaMe X 
belongs to the domain of attraction of an a-stable distribution if and only if 
the following two conditions hold: 

(i) the truncated moment function 

fi{x) := E{X'^1{\X\ ^ x}) (5.2) 
varies regularly with exponent 2 — a as x 00, i.e., 

fiix) - a;2-"Li(x), (5.3) 

where Li{x) varies slowly; 

(ii) either a = 2, or the tails of X are balanced: 

P(|X|>.) ^^+' -^00, (5.4) 

for some G [0, 1]. 

Proof Feller [3, Theorem XVII.5.2]. □ 

For the case a < 2, the following version is often more convenient. 

Theorem 5.3. Let < a < 2. A random variable X belongs to the do- 
main of attraction of an a-stable distribution if and only if the following two 
conditions hold: 

(i) the tail probability P{\X\ > x) varies regularly with exponent —a as 
X — >■ 00, i.e., 

P(|X| > x) - x-"L2{x), (5.5) 

where L2{x) varies slowly; 

(ii) the tails of X are balanced: 

F{X > x) 



¥{\X\ > x) 
for some p+ G [0, 1]. 



p+, x — >■ 00, (5-6) 



14 SVANTE JANSON 

Proof. Feller [3, Corollary XVII.5.2]. □ 

We turn to identifying the stable limit distributions in Theorems 5.2-5.3 
explicitly. 

5.1. The case a < 2. If the conditions of Theorem 5.2 or 5.3 hold for some 
a <2, then the conditions of the other hold too, and we have, by [3, (5.16)], 

2 — a 

L2{x) ^ LAx), X — >■ DO. (5-7) 

a 

Furthermore, by [3, (5.6)], with a„, 6„ as in (5.1) and M and A the canonical 
measure and Levy measure of the limit distribution £, 



l-OO 

n¥{X > Qnx) A{x, oo)= dM(y), x >0, (5.8) 

Jx 

and, by symmetry, 

/X 
y"^dM(y), x > 0. (5.9) 
-oo 

In particular, 

n¥{\X\ > an) ^ A{y : \y\ > 1} G (0, oo); (5.10) 

conversely, we may in (5.1) choose any sequence (a„) such that raP(|X| > a„) 
converges to a positive, finite limit. (Any two such sequences (a„) and (a^) 
must satisfy a„/a^ — >■ c for some c G (0, oo), as a consequence of (5.5).) 
If 

nP(|X| > a„) ^ C > (5.11) 

and (5.5)-(5.6) hold, then (5.8)-(5.9) hold with A(a;,oo) = p+Cx"" and 
A{—oo,—x) =p-Cx~°', where p_ := 1 —p+. Hence, (3.2)-(3.3) hold with 

c+=p_|_Ca, C-=p-Ca. (5-12) 

Consequently, the limit distribution is given by (3.6) where, by (3.11)-(3.12), 

7= (C7a(-r(-a)cos^))^/° = (Cr(l - a) cos ^)^/", (5.13) 
P=p+-p.. (5.14) 

For a = 1 we interpret (5.13) by continuity as 

7 = Cf. (5.15) 

Theorem 5.4. Let < a < 2. Suppose that (5.5) -(5.6) hold and that an 
are chosen such that (5.11) holds, for some C. Let 7 and /3 be defined by 
(5.13)^(5.14). 

(i) 7/0 < a < 1, then 

^ AS„(7,^,0). (5.16) 
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(ii) If 1 < a < 2, then 

>Sai7,/3,0). (5.17) 

On 

(iii) If a = 1, then 

^ILZl!^A^S,{j,(5,0), (5.18) 
where 7 is given by (5.15) and 

bn:=anEsm{X/an). (5.19) 

Proof. Feller [3, Theorem XVII. 5. 3] together with the calculations above. 

□ 

Example 5.5. Suppose that < a < 2 and that X is a random variable 

such that, as re — )• 00, 

P(X > x) ~ Cx-", (5.20) 

with 00, and F{X < -x) = o(a;-"). Then (5.5)-(5.6) hold with L2{x) := 
C and p+ = 1, and thus p_ := 1 — 1?+ = 0. We take a„ := n^/°; then (5.11) 
holds, and thus (3.2)-(3.3) hold with 

c+ = Ca, c_ = 0; (5.21) 

hence, (5.13)-(5.14) yield 

7 = (Cr(l - a) cos ^) (5.22) 

and /3 = 1. Consequently, Theorem 5.4 yields the following. 

(i) If < a < 1, then 

AS«(7,1,0). (5.23) 

The limit variable Y is positive and has by Theorem 3.12 and (5.22) 
the Laplace transform 

Ee"*'*' =exp(-CT(l-a)r), Ret ^ 0. (5.24) 

(ii) If 1 < a < 2, then 

^S„(7,l,0). (5.25) 

The limit variable Y has by Theorem 3.12 and (5.22) the finite Laplace 
transform 

Ee-*^ = exp(C|r(l-a)|t"), Ret^O. (5.26) 

By (4.6) and (5.22), the density function fy of the limit variable sat- 
isfies 

/(O) = C-^/°^|r(l - a)r^/'*|r(-l/a)r^ (5.27) 
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(iii) If a = 1, then 

^?^i^ = ^-6„ ASi(7,l,0), (5.28) 
n n 

where, by (5.15), 7 = C7r/2 and 

bn:=nEsm{X/n). (5.29) 
We return to the evaluation of 6„ in Section 5.2. 

Example 5.6. Suppose that < a < 2 and that X ^ is an integer-valued 
random variable such that, as n — )■ 00, 

P(X = n) - cra-"-^ (5.30) 

Then (5.20) holds with 

C = c/a (5.31) 
and the results of Example 5.5 hold, with this C. In particular, (5.22) yields 

7" = -cr(-a)cos^, (5.32) 
and both (5.24) and (5.26) can be written 

Ee-*^ = exp(cr(-Q)t"), Ret ^ 0; (5.33) 

note that T{-a) < for < a < 1 but T{-a) > for 1 < a < 2. 
Taking t imaginary in (5.33), we find the characteristic function 

Ee'*^ = exp(cr(-a)(-ii)") = exp(cr(-a)e-'^g^(*)'^"/2|tr) , t G M. 

(5.34) 

5.2. The special case a = 1. Suppose that, as x — >■ 00, 

¥{X > x) ~ Cx-'^ (5.35) 

and F{X < —x) = o{x~^), with C > 0. Then Example 5.5 applies, and 
(5.28)-(5.29) hold. We calculate the normalising quantity 6„ in (5.28) for 
some examples. 

Example 5.7. Let X := 1/U, where U ^ U(0, 1) has a uniform distribution. 
Then F{X > x) = x'^ for x ^ 1 so (5.35) holds with C = 1 and (5.15) yields 
7 = 7r/2. Furthermore, X has a Pareto distribution with the density 

10, X ^ 1. 

Consequently, by (5.29), 

roo poo 

bn = nsm{X/n) = n / sin(a;/n)x~^ dx = / sm{y)y~^ dy 

Jl Jl/n 

siny — y f°° siny 

= log n + / 2 dy + / — 2- 

Ji/n y Jl y 

/■■^ sin?/ - ■yljy < 1| , 

= logn+ / ^ — ^dy + o 1 =logn + l-7 + o 1 , 

Jo y 
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where 7 is Euler's gamma. (For the standard evaluation of the last integral, 
see e.g. [7].) Hence, (5.28) yields 

n 

or 

n 



:^_(logn + l-7) ASi(7r/2,l,0). (5.37) 
^-logn ASi(7r/2,l,l-7). (5.38) 



Example 5.8. Let X := 1/Y, where Y ~ Exp(l) has an exponential dis- 
tribution. Then F{X > x) = 1 — cxp(— 1/.t) ~ x^^ as a; — > oo so C = 1 
and (5.15) yields 7 = 7r/2. In this case we do not calculate bn directly from 
(5.29). Instead we define U := 1 — and X' := 1/U and note that U has 
a uniform distribution on [0, 1] as in Example 5.7; furthermore 

1 1 - l+Y 

This is a positive random variable with finite expectation 



e V 



E{X'-X)= —^e-ydy= I- dy = 7, 

Jo [l-e y}y Jo \l-e y y J 

(5.40) 

see e.g. [2, 5.9.18] or [7]. 

Taking i.i.d. pairs {Xi,X'^ = {X,X') we thus have, with S'^ := ELi^i' 
by the law of large numbers, 

E{X' - X) = 7. (5.41) 



S'n — Sn p 



n 

Since Example 5.7 shows that S'^/n — log ri Si(7r/2, 1, 1 — 7), it follows 
that 

5„/n-logn ASi(7r/2,l,l-27). (5.42) 

We thus have (5.28) with 

6„ = logn + 1-27 + 0(1). (5.43) 

5.3. The case a = 2. If a = 2, then a„ in (5.1) have to be chosen such 
that 

n/x(a„) 



i2 



C (5.44) 

for some C > 0, see [3, (5.23)]; conversely any such sequence (a^) will do. 

Theorem 5.9. // fj,{x) is slowly varying with fj,{x) 00 as x 00 and 
(5.44) holds, then 

^ILZl^ N{0,C). (5.45) 

Proof. Feller [3, Theorem XVII.5.3]. □ 
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Example 5.10. Suppose that a = 2 and that X is a random variable such 
that, as a; — >■ oo, 

P(X > x) - Cx-^, (5.46) 
with C > 0, and P(X < -x) = o(x"2)^ -phcn (5.4) holds with p+ = 1, and 
thus p- := 1 — p+ = 0. Furthermore, as x ^ oo, 

li{x)=¥.{^j 2tdtl{\X\^x}^ j l{t^\X\^x}2tdt 

PX PX 

= 2tF{t ^ \X\ ^ x) dt = 2tF{\X\>t)dt-x^F{\X\> x) 
Jo Jo 

= (l + o(l))^ 2iCt-^dt + 0(l) ~ 2Clogx. (5.47) 

Thus (5.3) holds with Li{x) = 2Clogx. 

We take a„ := Vnlogn. Then //(a„) ~ 2C5logn = Clogn, so (5.44) 
holds and Theorem 5.9 yields 

^zl£^A^NiO,C). (5.48) 
V n log n 

6. Attraction and characteristic functions 

We study the relation between the attraction property (5.1) and the char- 
acteristic function (px{t) of X. For simplicity, we consider only the common 
case when a„ = n^/". Moreover, for simplicity we state results for (px{t), 
t > only, recalling (3.15) and ^px{0) = 1- 

Theorem 6.1. Let < a ^ 2. The following are equivalent. 

(i) ^ Z for some non- degenerate random variable Z. 

(ii) The characteristic function ipx of X satisfies 

cp^(t) = 1- {K-iX)t" + o{t'^) ast\0, (6.1) 

for some real k > and X. In this case, Z is strictly a-stahle and has 
the characteristic function (3.16). (Hence, |A| ^ Ktan^.J 

Proof. If (i) holds, then for every integer m, 

g ^ in ^ n ^ rn 

= ^T. ^Y.^ik-i)n+j ^^Y.^k. as n ^ oo, 

^ ' k=\ 3=\ k=\ 

with Zfc = Z i.i.d. Since also {mn)~^^°'Smn we have m~^/" X^feLi = 

Z, and thus Z is strictly ct-stable. 

We use Corollary 3.8 and suppose that Z has characteristic function 
(3.16). Then the continuity theorem yields 

<^x(t/n^/T^¥'z(t) = exp(-(AC-iA)r), i ^ 0; (6.2) 
moreover, this holds uniformly for, e.g., ^ t ^ 1. 
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In some neighbourhood {—to, to) of 0, ipx ^ and thus (px{t) = e'^^^^ for 
some continuous function ip : {—to, to) — >■ C with ■0(0) = 0. Hence, (6.2) 
yields (for n> 1/to) 

exp (ntp ( ^ . ) + (k — iA)t") = 1 + o(l), as n ^ oo, 



uniformly for ^ t ^ 1, which implies 



nil) ( ^ , \ + (k — \\)t°' = o(l), as n — >■ oo, 



since the left-hand side is continuous and for f = 0, and thus 
( t \ 

V' I , + (k — iA) — = o(l/n), as n ^ oo, (6.3) 
yn'-l'^ J n 

uniformly for ^ t < 1. 

For s > 0, define n := [s~°J and t := sTn}/"^ G (0, 1]. As s \ 0, we have 
n — >■ oo and (6.3) yields 

V,(s) = -{k - i\)s'^ + o(l/n) = -{k - i\)s'^ + o(s°). (6.4) 

Consequently, as s \ 0, 

^x{s) = e'^(^) = 1-{k- iA)s° + o{s'^), (6.5) 

so (6.1) holds. 

Conversely, if (6.1) holds, then, for t > 0, 

E e'*^"/"'^" = ifx (i/n^/")" = (l - (k - iA + o(l)) ^) " ^ exp (-(k - iA)^) , 

as 77 — )■ oo, and thus by the continuity theorem Sn/n^^'^ — >■ Z, where Z has 
the characteristic function (3.16). □ 

For a = 1, it is not always possible to reduce to the case when 6n = in 
(5.1) and the limit is strictly stable. The most common case is covered by 
the following theorem. 

Theorem 6.2. The following are equivalent, for any real b. 

1^ d 

(i) — — 61ogn — > Z for some non-degenerate random variable Z. 
n 

(ii) The characteristic function (px of X satisfies 

ipx {t) = l-{K- iX)t - ibt log t + o{t) ast\0, (6.6) 

for some real k > and A. In this case, Z is 1-stable and has the 
characteristic function (3.23). (Hence, \b\ ^ 2k/tt.) 

Proof, (ii) (i). If (6.6) holds, for any k G M, then, as t\0, 

log ipx{t) = -{K-iX + o{l))t - ibtlogt (6.7) 
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and thus, as n — > oo, for every fixed t > 0, 

= expfnf— (k — iA + oil))— — ib— log — ) — i6t logn^ 
\ \ n n nJ J 

exp (— (k — \\)t — ibt log t) 

which shows (i), where Z has the characteristic function (3.23). 

Furthermore, for use below, note that (3.23) implies |(^z(i)| = e"*^* for 
t > 0. Since |99^(t)| ^ 1, this shows that k ^ 0. Moreover, if k = 0, then 
|¥'z(i)| = 1 for i > 0, and thus for all t, which implies that Z = c a.s. for 
some c G M, so Z is degenerate and 6 = 0. Hence, (6.6) implies k ^ 0, and 
K = is possible only when 6 = and Sn/n A. 

(i) =^ (ii). Let 71 := |6|7r/2 and /3i := — sgn6. Let Y and Yi be i.i.d., 
and independent of (Xj)^ and Z, with distribution 81(71, 0). (If 6 = 
we simply take Yi := 0.) Then Yi has, by (3.6), the characteristic function 

(/?y(i) = exp(-7it + i6ilogi), t > 0. (6.8) 

By Theorem 3.3(ii), 

n 

J2Yi = nY-bn log n. (6.9) 

i=l 

Define Xi := Xj + Yi. Then, 

iVx, = iVx, + iVyi^^ + r-61ogn Az + y. (6.10) 
n ^-^ n ^-^ n ^-^ n 

1=1 i=l 1=1 

Thus, by Theorem 6.1, for some ^2 > and A2, 

(pxit)(pY{t) =Ee'*^' = 1- {K2-iX2)t + o{t) asi\0, (6.11) 
and hence, using (6.8), 

ipx{t) = Ee^*^7</?y(t) = 1 - (k2 - iA2 - 7i)i - ibtlogt + o{t), (6.12) 

which shows (6.6), with k = K2 — 71 € M. 

Finally, we have shown in the first part of the proof that (6.6) implies 
K> 0, because Z is non-degenerate. □ 

We can use these theorems to show the following. 
Theorem 6.3. Let < a ^ 2. Suppose that X is such that 

n 

-1/a -.^ d . 



n 

1=1 



Y^Xi^Z, (6.13) 
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where Z is an a-stable random variable with characteristic function (3.16) 
and that Y is a random variable withEY"' < oo. Let be indepen- 

dent copies ofY that are independent of {Xi)^ . Then 

n 

J2 XiYi A Z' := (E y°) ^^"Z, (6.14) 

i=l 

where the limit Z' has the characteristic function 

(^z'(*) = exp(-(Ey"K-iEF"A)t"), t ^ 0. (6.15) 

// Z~ S„(7,/3,0) (where ^ = if a = 1), then ^' ~ Sa((EF")V«7, /3, 0). 

Proof. By Theorem 6.1, for t ^ 0, 

¥Jx(i) = l-(^c-iA)r + rr(t), (6.16) 

where r{t) as t \ 0. Furthermore, (6.16) imphes that r{t) = 0{1) as 
t -)• oo, and thus r{t) = 0(1) for t ^ 0. 

Consequently, for t > 0, assuming as we may that Y is independent of X, 

ipxYit) = Ee'*^^ =Eipx(tY) = E(l - (k - iA)ry" + t'^Y'^r{tY)) 

= 1- (/«-iA)rEy" + rE(y"r(iy)), (6.17) 

where E[Y°'r{tY)) as t \ by dominated convergence; hence 

<^^y(t) = 1 - (K-iA)t"Ey" + o(r) ast\0. (6.18) 

Theorem 6.1 appHes and shows that n^^/" "121=1 ^i^i ~^ ^' ^ where Z' has 
the characteristic function (6.15). Moreover, by (3.16), (Ey")^/" has this 
characteristic function, so we may take Z' := (Ey")^/". 

The final claim follows by Remark 3.6. □ 

Theorem 6.4. Suppose that X is such that, for some real b, 

n 

n-^^Xi-blogn^ Z, (6.19) 

i=l 

where Z is a 1-stable random variable, and that Y ^ is a random vari- 
able with Eyiogy < oo. Let (ii)i° be independent copies of Y that are 
independent of {Xi)f. Then, with /j, :=EY, 

n 

n-^ J2 XiYi - bn log n A Z' := //Z - 6(E(y log Y) - log fi) . (6.20) 

i=l 

Z has the characteristic function (3.23) for some k and X, and then the 
limit Z' has the characteristic function, with v := E(y logy), 

ipz'{t) = exp{-{nK-\-i{bu - n\)t) -ibntlogt), t > 0. (6.21) 

Si(7,/3,(5), then Z' ~ Si{n^ , P , fiS - bu) . 
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Proof. By Theorem 6.2, for t ^ 0, 

ipx{t) = 1 - (k- iA)t - ibfAogt + tr{t), (6.22) 

where r{t) — )■ as t \ 0; moreover Z has the characteristic function (3.23). 
Furthermore, (6.22) impUes that r{t) = O(logt) as t ^ oo, and thus r{t) = 
0(1 + log_^ t) for t ^ 0. 

Consequently, for t > 0, assuming as we may that Y is independent of X, 

^XY{t) = E^x{tY) 

= 1 - {K-[X)tKY - ibtK{Ylog{tY)) +tE{Yr{tY)), 
= 1 - {fiK-ifiX + ibE{YlogY))t-ibfitlogt + tE{Yr{tY)), 
where M(Yr{tY)^ — t- as t \ by dominated convergence; hence 

fxrit) = 1 — {P'K — ifiX + ibu)t — i6//ilog t + o{t) as t \ 0. (6.23) 

Theorem 6.2 apphcs and shows that X^ILi ^lYi — bji log n — ^ Z' , where 
Z' has the characteristic function (6.21). Moreover, it follows easily from 
(3.23) that fiZ — b(E,{Y log Y) — jj. log ^) has this characteristic function, and 
thus (6.20) fohows. 

Finally, if Z ~ 81(7, /J, 5), then b = -/3^ by Remark 3.10 and it follows 
easily from Remark 3.6 that Z' ~ Si(//7, P, /id — bu); alternatively, it follows 
directly from (6.20) and (3.6) that Z' has the characteristic function 

(pz'{t) = <fzilJ't)exp{-ibt{i' - /xlog/x)) 

= exp (^-7At|i I (1 + i^- sgn(t) log |i |^ + id fit - ibtu^ . (6.24) 

□ 

Example 6.5. Let X := U/U', where U, U' U(0, 1) are independent. By 
Example 5.7 and Theorem 6.4, with Z ^ Si(7r/2, 1, 1 — 7), b = 1, /j, := EU = 
1/2 and 

/•I 1 

1/ := EC/logC/ = / xlogxda; = — -, (6.25) 
Jo 4 

we obtain 

logn — > -Z — V -\ — log - = -Z -\ log 2 Si — , 1, 

n 2^ 2 2^2 2 42^ V4' '4 2 

(6.26) 

Example 6.6. Let X := Y/Y' where Y,Y' ~ Exp(l) arc independent. 
(Thus X has the i^-distribution i*2,2-) By Example 5.8 and Theorem 6.4, 
with Z ~ Si(7r/2, 1, 1 - 27), 6 = 1, := EF = 1 and 



v-EY\ogY=l xloga;e"^dx = r'(2) = 1-7, (6.27) 
Jo 



we obtain 



— -logn^Z-i/ = Z- l + 7~ Si(7r/2, 1, -7). (6.28) 

n 
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This is in accordance with Example 5.7, since, as is weh- known, U := 
Y'/{Y+Y') ~ U(0,1), and thus we can write X = {¥+¥')/¥'-! = 1/U-l. 

Example 6.7. Let X := V^/W where V ~ U(-i, i) and W ~ Exp(l) are 
independent. By Example 5.8 and Theorem 6.4, with Z ~ Si(7r/2, 1, 1 — 27), 
b=l, fj,:=EV'^ = 1/12 and 

fi/2 r^3 ^31 1/2 

1/ := 2Ey2log|F| = 4 / a;^logxda; = 4 
Jo 

31og2 + 1 



— log X 

3 ^ 9 



18 



we obtain 

Sn 1 



^1 d 1„ ^1, 1 /TT 5-67 + 61og2 

--^logn^-Z-.+-log-~Si(-,l 



24' 



36 



(6.29) 
). (6.30) 



Equivalently, using Remark 3.6, 
24S'„ 2 , d 2 „ 24i/ 2 



- - logn A -Z log 12 - Si 1, 1, - - - 27 + log - 

TT TTTTTT \7r\3 



6- 

(6.31) 

This is shown directly in Heinrich, Pukelsheim and Schwingenschlogl [5, 
Theorem 5.2 and its proof]. 

Example 6.8. More generally, let X := V"^ /W where V ~ \J{q — 1, 1) and 
W ~ Exp(l) are independent, for some fixed real q. By Example 5.8 and 
Theorem 6.4, with Z ~ Si(7r/2, 1, 1 - 27), 5 = 1, 

2 / 1n2 1 3g2_3g_^l 



^ = EF2=(^Ef) +Vary=(g--) +- = 



(6.32) 



and 



1/ := 2EF^log|F| = 2 / x^log|x|dx = 2 



^ x'^ X 
ylogN-- 



3n9 



log |g| + (1 - qf log |1 - q\ ^ 3g^ - 3g + 1 



we obtain 

n 

Equivalently, using Remark 3.6 



/i log n /xZ - 1/ + ^ log ' 1' (1 ~ 2')')^ ~ 



(6.33) 

(6.34) 



^ log n — > Z h log /i — -^^ 



with 

L 2 log i^^i + (1 - (^)'^ log |1 - q\ 
^-=3-'^-' 3g^-3g + l + 3 



^ Si (^,1,6,), (6.35) 

3rr' - 'M + 1 1 



12;u" 
(6.36) 

This is shown (in the case ^ g ^ 1) directly in Heinrich, Pukelsheim 
and Schwingenschlogl [6, Theorem 4.2 and its proof]. 
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